Spin-chain description of fractional quantum Hall states in the Jain series 
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We discuss the relationship between fractional quantum Hall (FQH) states at filling factor v = 
p/(2p+ 1) and quantum spin chains. This series corresponds to the Jain series v — p/(2mp+ 1) with 
m = 1 where the composite fermion picture is realized. We show that the FQH states with toroidal 
boundary conditions beyond the thin-torus limit can be mapped to effective quantum spin S = 1 
chains with p spins in each unit cell. We calculate energy gaps and the correlation functions for both 
the FQH systems and the corresponding effective spin chains, using exact diagonalization and the 
infinite time-evolving block decimation (iTEBD) algorithm. We confirm that the mass gaps of these 
effective spin chains are decreased as p is increased which is similar to S = p integer Heisenberg 
chains. These results shed new light on a link between the hierarchy of FQH states and the Haldane 
conjecture for quantum spin chains. 

PACS numbers: 71.10.Pm, 75.10.Kt, 73.43.Cd 



I. INTRODUCTION 

The fractional quantum Hall (FQH) stated interact- 
ing cold electrons in two-dimensional space in a strong 
perpendicular magnetic field, exhibits fascinating phases 
with fractionalized excitations and topological order^r— 
Ever since its discovery three decades ago, the FQH sys- 
tem has inspired a huge amount of experimental and the- 
oretical effort, due to its richness in phenomenology and 
mathematical structure. New developments include the 
observation of the FQH in graphenei, topological quan- 
tum computing 8 , and systems of rapidly rotating bosons 
which arc formally very similar to those of an electron 
gas in a magnetic field. 9 

On the other hand, it has been pointed out that the 
hierarchy of FQH states has striking similarities to the 
quantum spin chains. 10 Haldane conjectured^ 1 - that half- 
integer SU(2) Heisenberg chains support gapless excita- 
tions, protected by a topological term in the effective 
action, while the integer spin chains develop a mass gap. 
A similar structure appears in the FQH effect. At filling 
factors v < 1, quantized conductance plateaus only occur 
at rational v with odd denominator, while in the vicinity 
of euen-denominator fractions metallic behavior is sus- 
tained. Hence, it is important to establish whether the 
similarities arc merely accidental or whether the struc- 
ture of low-energy excitations in these systems has a re- 
lated microscopic origin. 

Recently, a framework for studying this connection 
was proposed which maps FQH systems with torus ge- 
ometry to one-dimensional (ID) discretized models. It 
was realized that universal features of many FQH phases 
are retained in its thin-torus (or Tao-Thouless^ TT) 
limit ; 13 ' 14 where the interacting system becomes a triv- 
ial ID charge-density-wave (CDW) state. FQH states at 
odd-denominator filling fraction can be deformed into the 
TT limit without closing the energy gap, as has been rig- 
orously shown at the Laughlin fractions l/qM-~— Based 
on this property, the v = 1/3 FQH state on a torus be- 
yond the TT limit has been mapped to an 5 = 1 spin 



chain, and it is shown that the ground state is a gapful 
state which is adiabatically connected both from the Hal- 
dane gap phase and the large- D phased - — This special 
situation is realized due to the breaking of the discrete 
symmetries of the effective spin model. On the other 
hand, notably different behavior is found in states at 
even-denominator filling. For example, a gapless state 
at filling fraction v — 1/2 undergoes a phase transition 
from a gapped TT state to a gapless phase upon deforma- 
tion of the torus. This can be interpreted by mapping the 
system to the 5 = 1/2 XXZ spin chain which undergoes 
a phase transition from the ferromagnetic state to the 
Tomonaga-Luttinger liquid phasei 13 ' 16 From the above 
results, we speculate that FQH states with odd(even) 
denominator filling fractions are related to integer(half- 
integer) 5 spin chains. 

In this article, we extend this approach to more gen- 
eral FQH states. As discussed by Jain, 5 FQH states 
at v = p/(2mp + 1) can be described as the compos- 
ite fermion picture, where a state in which 2m quantum 
fluxes are attached to noninteracting electrons of the pth 
Landau level is projected onto the lowest Landau level. 
Among these Jain series, we turn our attention to m = 1 
cases v — p/(2p + 1), since this series is very impor- 
tant to connect v = 1/3 (p = 1) and v = 1/2 (p — > oo) 
systems, and also straightforward extensions of the above 
spin mapping are possible. 18 Therefore, we consider map- 
ping of these states to quantum spin chains and study 
their properties. 

The rest of this paper is organized as follows. In Sec.HIl 
we explain how the FQH states with torus geometry are 
described by ID discretized models. In Sec. IIIIl we inves- 
tigate how the v =p/ (2p+ 1) FQH states are mapped to 
spin variables, and conclude that the effective Hamiltoni- 
ans are 5=1 quantum spin chains with p-site unit cells. 
In Sec. IIV| we analyze the properties of effective spin 
chains numerically, using exact diagonalization and infi- 
nite time-evolving block decimation (iTEBD) algorithm, 
and show that the effective model for v = p/(2p + 1) 
behaves like an 5 = p quantum spin chain. Concluding 
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remarks are given in Sec. [Vj In appendices, we present 
mathematical proofs for lemmas used in the spin map- 
ping. 



II. ID DESCRIPTION OF THE FQH STATES 

We consider a model of N interacting electrons on a 
torus with circumference Li in xi direction (I = 1,2). 
When the torus is pierced by N s magnetic flux quanta, 
L\Li — 2irN s is s atisfied , where we have set the mag- 
netic length Ib = y/H/eB to unity. In the Landau gauge, 
A = BX2X1, a complete basis of N s degenerate single- 
particle states in the lowest Landau level, labeled by 
k = 0, . . . , N s — 1, can be chosen as 



ipk 



1 



j(k 1 +nL 2 )-^(x2+k 1 +nL2) 2 



(1) 



where k\ = 2nk/Li is the momentum along the in- 
direction. In this basis, any translation-invariant two- 
dimensional Hamiltonian with two-body interaction as- 
sumes the following ID lattice model: 



H = V km ; 

\m\<k<N B /2 



Vkr, 



(2) 



+k C j+k+m C ] 



where the matrix element V Km specifies the amplitude of 
a pair-hopping process. In this model, two particles sep- 
arated k + m sites hop to steps to opposite directions, 
then their distance becomes k — m sites (note that m 
can be or negative). The to = terms can be re- 
garded as the electrostatic repulsion. At filling fraction 
v = p/q < 1) the Hamiltonian commutes with the center- 
of-mass magnetic translations, T}, along the cycles. They 
obey T x f 2 = e 27rip /«f 2 Ti, so that the operators H, T\, 
and T 2 q commute each other. From the periodic bound- 
ary conditions, TJ " = 1, two conservation numbers are 
given as 



Ti 



,27riKi 



rpq 

1 2 



(3) 



All energy eigenstates are (at least) g-fold degenerate, 
and all states can be characterized by a two-dimensional 
vector Ki = 0, 1, • • ■ ,N s /q — 1. K± denotes center-of- 
mass quantum numbers for the x 2 direction. 

For small L\ the overlap between different single- 
particle wave functions ([T]) decreases rapidly and the ma- 
trix elements V km are simplified considerably. As L\ — >• 
one finds that 



Vk; 



o 2 2 / j- 2 

V k0 e- 2 * m ' L i 



(4) 



thus the to ^ terms are exponentially suppressed for 
generic interaction in this limit. The remaining (m = 0) 
problem becomes trivial: Ground states at any v = p/q 



are gapped periodic crystals (with a unit cell of p elec- 
trons on q sites) and the fractionally charged excita- 
tions appear as domain walls between degenerate ground 
states. This is the state that Tao and Thouless proposed 
to explain the quantum Hall effect^ therefore we of- 
ten refer to this limit as the Tao-Thouless, or thin-torus, 
limit. In the TT limit the Hamiltonian can be written as 



Utt = /J Vko 



nj+ k nj, 



(5) 



where n,- 



d fa- 



The ground state of Eq. §5§ is ap- 
parently a CDW state with g-fold degeneracy, since the 
electrons favor being located as far as possible from each 
other. In addition we can interpolate between the solv- 
able limit and the bulk by continuously varying a single 
variable, L\. 

In this paper we consider a truncated Hamiltonian of 
([2]) with only the two most dominant electrostatic terms 
and one hopping term, 



H it = } y [V w hjn j+ i + V 2 oh 3 h j+2 



j 



+ ^21(4+14+2^+3^ + H.C.)] . 



(6) 



This provides a good approximation of a short-range in- 
teraction. We consider a Trugman-Kivelson type pseu- 
dopotential V(r — r') cx V 2 i5(r — r')^ on a thin torus 
(Vio > 2V20 3* others), where the matrix elements for 
L-2 — > 00 are 



V km oc (k 2 - m 2 )e- 2 * 2 ( k2+m2 ^ L l 



(7) 



For Coulomb interaction, the longer range electrostatic 
terms V k o are nonnegligible. 



III. SPIN MAPPING OF m = 1 JAIN SERIES 

In order to study properties of the FQH states de- 
scribed by the ID model ([S]), we consider spin mapping 
of this system. This mapping is also interesting to know 
the relationship between different physical systems: FQH 
states and quantum spin chains. 

In the case v = 1/2 {p — >-oo,a non-FQH state) and 
v = 1/3 (p = 1), one can find trivial subspaces of the 
truncated Hamiltonian ([5]) for the spin mapping. For 
v = 1/2 , a subspace of the full Hilbert space is required 
where each pair of sites (2n — l,2n) has one particle. 
Therefore, there are only two possible states for a pair of 
sites in the subspace, and they can be related to S = 1/2 
spin variables as |01) -> |f) and 1 10) -> |4_). Thus the 
system is mapped to an S = 1/2 spin chain^ 

For v = 1/3, the TT limit of the system 
is the threefold-degenerate CDW state I^tt) = 
|- • ■ 010 010 ■ • •), where the underlines denote unit cells. 
The degenerate states have different center-of-mass quan- 
tum numbers {K\) and there are no matrix elements 
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between them. In this case, the subspace is given by 
configurations generated by applying t/ 21 several times 
to |^tt)- Then the spin variables can be introduced as 
|010) -> |o), |001) -> |+>, 1 100) -> |-). This 5 = 1 model 
has been discussed in previous work>i^ ? i2£ 

In this paper, we consider the extension of the above 
spin mapping. As discussed by Jain^S FQH states at 
v = p/(2mp + 1) can be described as the composite 
fermion picture. For the filling fractions with m = 1 
Jain sequences, v — p/(2p + 1), the CDW state in the 
TT limit is given by | - • • 0(10) p 0(10) p ■ • •) where the unit 
cell 0(10) p denotes a configuration which consists of 
and p times of 10. This clearly minimizes repulsion of 
the electrostatic terms. Since these states have one or 
two particles in every three sites, a natural extention 
of the spin mapping of the truncated model ([6]) to the 
v = 1/2,1/3 FQH states is expected. In what follows, 
we discuss how the subspace of the m = 1 Jain series is 
identified in term of spin variables. 



Here, we choose a unit cell for the ground-state wave 
function in the TT limit as |- ■■ 01010 01010 ■■ ■). This 
state has two particles in each five-site unit cell. Us- 
ing the property of Eq. ([lip we can immediately confirm 
that the states ■ ■ ■ 000?? ■ ■ ■) and |- ■ ■ ??000 ■ ■ ■) where 
? = 1,0 do not appear in the subspace. The states 
|- •■ 001001 ■••) and |- •• 1 00100 •• •) which can generate 
the states ■■ - 00011 •■•) and ■ ■ ■ 11000 - ■ ■) are not in- 
cluded in the subspace either. Since the vanishing states 
| - • • 00100 • • •} can only be generated from the states 
|- •• 00100 1 •••} or |--- 100100- •■), they cannot be ele- 
ments of the subspace either. Therefore each unit cell 
always includes two particles, so that the states are iden- 
tified as 5 = 1 spin variables by inserting appropriately 
between the two l's such as 



|01010) -> |01[00]10) -> |oo) , (12) 
|00110) -> |00[10]10) -> |+o) , (13) 
|01100) -> |01[01]00) -> |o-) , (14) 



A. The subspace 

Let us now consider extensions of the above 5 = 1 
mapping to other Jain series. We define a local operator 
which gives a pair hopping process in V21, 



Uj - C ] + l C ]+2 C J+3 C j- 



(8) 



Moreover, we introduce a state vector |^ r ) where all 
electrons do not occupy the nearest two sites (e.g., 
• • • 01010010 •••)). Then it satisfies the following rela- 
tion: 



U] |* r > = 0. 



(9) 



We choose such a state vector as a "root state" of general 
configurations. In this case, all configurations in that 
subspace can be written as 



*«= b uuh,~] l%,j»,-]> ; 

\*uuh,---])=nu jh \v r )=p {juj2t ... ] \* r ), 



(10) 



where by lt j 2t ,,.-\ are expansion coefficients for the config- 
urations of the state \^[j 1 ,j 2 ,---]) ( see Appendix A 
process which generates a configuration \^[j 1 ,j 2 ----]) fr° m 
the root state |\Ev) is specified by the series [ji, j'2, ■ • •]■ 
Therefore, for each configuration of states, their creation 
operators P[j lj2j ...] are defined. Using the basis of (ITU1) . 
we find the following condition around a site j, 



\* r ) = |. - -010 • - •) 



I* 



[j'l ,32 



.]) ^ |- - -000- - •) , (11) 



where the site j is the center of the underlined three sites. 
Proof of Eq. pT|) will be given below. 

Before presenting results for general v = p/(2p + 1), 
let us consider a case of v = 2/5 as the simplest example. 



and so on. Thus the subspace of the truncated Hamilto- 
nian for v = 2/5 can be mapped to two 5 = 1 spin 
variables just like the case of v = 1/3. 

In fact, the property ([TTj) is a special case (p = 1) of 
the lemma that the number of electrons in a unit cell 
should be unchanged: 



|tf r ) = |--- 0(10) p , ---) 

=> J2 (*[j I j 2 ,---]N*tfi, J v--]> >p«. ( 15 ) 
1=3 

where site j is the first site of the underlined 2p s + 1 
sites, and the root state |\fv) does not include | 111) (see 
Appendix [Bj) . 

At v — p/(2p + 1), wc choose the root state \^ r ) 
as | - • • 0(10) p 0(10) p • • •). Then using the condition (fT5|i 
with p s = p and the particle conservation, the local par- 
ticle conservation in each unit cell can be confirmed. We 
now can use the condition (|15p again with p s = p — 1 to 
confirm that there is only one electron in sites j and j + l. 
By performing this process from p s = p to p s — 1 in simi- 
lar ways, we can confirm that the truncated Hamiltonian 
© can be mapped to 5 = 1 quantum spin chains (see 
Fig. [T]) by defining spin variables like the case v = 2/5. 

We can introduce more general states which are not 
generated from the simple root state in Eq. dl~5l . but 
from more complicated configurations. However, these 
states can also be decomposed into domains with unit 
cells of v = p/(2p + 1) with different p. For example 
• • • 01001010010 • • •) can be given by p = 1 and p = 2 
unit cells |- - - 010 01010 010 ■■ ■). We expect that those 
states have higher energy than the states without do- 
mains, because they have larger unit cells. 
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- (5 = o H q-Q )- -( o-o-c^><><)) - 

y = 1/3, 2/3 ^ = 2/5,3/5 v = 3/7, 4/7 



FIG. 1. Effective S = 1 spin chains for the v = p/(2p + 1) 
FQH states. There are p spins in each unit cell. 
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TABLE I. Relationship between fermion and 5=1 spin op- 
erators for the v = 2/5 FQH state. 



B. The effective Hamiltonian 



neighboring unit cells, the relations are 



In the case v = 1/2, as discussed in Ref. \L3L the trun- 
cated Hamiltonian ([6]) exactly gives an 5 = 1/2 XXZ spin 
chain. In the subspace of this system, S = 1/2 spin oper- 
ators are introduced as c\ n c 2n+1 -)• S~, c\ n+1 c 2n £+, 



1/2 + 5*, a 



Eq. ([5]) becomes 



2n+l u 2n+l 



1/2 



Thcrfore 



%XXZ — ^(^n 'S'n+l + ^n+l) + 



where A = (2V20 — Vio)/4V2i. This model undergoes a 
phase transition from a ferromagnetic phase to a gaplcss 
XY phase at L\ ~ 5.3. 

For ^ = 1/3, in addition to the mapping in the 
preceeding worksji^— we have to introduce a function 
f s (S* +1 — 5 T Z j) which stems from the electrostatic terms 



f s (x) = V 3 - x ,o, 



(17) 



where Vko — (k > 2) in the truncated Hamiltonian ((6|). 
The amplitude of electrostatic terms depends only on the 
difference between the S z of neighboring two spins. Then 
the effective spin Hamiltonian is obtained as 



+ /s('S'n+l ^n) 



(18) 



In the cases of v = p/ (2p + 1) with p > 2, we consider a 
unit cell which consists of p spins as shown in Fig. [TJ Wc 
introduce spin operators of the jth site in the nth unit cell 
as S" n . We should note that the spin-spin interactions 
inside of a unit cell are different from those involving the 
neighboring two unit cells. We consider the roots with 
A'i = which give the ground states. For the hopping 
process (V21) in a unit cell, the corresponding S = 1 spin 
operators are related as 

C(2p+l)^n(2p+l) + l ~* 2 1/2 ^l,n*Sl,n' 0-9) 
^!(2p+l)+2p-l^ri(2p+l)+2p ~ 2 ^ Sp,nSp,ni (20) 

and their Hermitian conjugates. For those between 



C(2p+l)+2jSi(2p+l)+2j-l 



_v _ 0-3/2 nz o+ c+ C— 

°j-l,n 7-l,n a 7,n a 7",n' 



C n(2p+l)4-2j C n(2p+l)+2j + l 



— 3/2 OS Q — Q— C+ 



(21) 



(22) 



We also need to introduce contributions from the electro- 
static terms between unit cells as a function fi(S^ +1 — S*) 
defined by 



fi(x) = V 2 - X: o- 



(23) 



As a simplest example, we consider the v = 2/5 state 
again. In this case there are two spins in each unit cell 
(S",T"). The relation between combinations of original 
fermion operators and the corresponding 5 = 1 spin op- 
erations are summarized in Table|H Then the spin Hamil- 
tonian can be written as 



E 



21 /rp-\-rpZ QZ q- j_ JiZji- o+ QZ 

2 V n n^n^n ' n n °n 

-S^S+t +1 f- +1 + R.c.) (24) 



The effective Hamiltonian of the m = 1 Jain states 
with general p is given by 




V°i,n°j,n'-'j+l,n' ; 'j + l, 



J," j,n ^j+l,n Sj+l,n) + H.C. + fi(Sj + ln Sj >n ) 



(25) 



The effect of the electrostatic terms can be neglected 
if the functions in the last terms are approximated as 
f s ,i(x) oc x, since the spin variables cancel each other 
and these terms only give constants. 



IV. NUMERICAL ANALYSIS 

In order to confirm the validity of our spin mapping of 
the FQH state, we numerically calculate energy gaps of 
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FIG. 2. Energy gaps of the original FQH systems (a), and those of the corresponding effective spin chains (b). L is the system 



size which describes the number of electrons in FQH systems (L/u = 15 ~ 27, 10 ~ 25, 14 

~ 15). 
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FIG. 3. Correlation functions of the effective spin chains with (a) and without (b) static terms obtained by iTEBD, while 
(a') and (b') are their logarithms. S z (n) is the total spin of the n-th unit cell S z (n) — X^j™np+i &j,n- r denotes the distance 
between two total spins where the lengths of unit cells are set to be unity. 



v = 1/3,2/5,3/7 FQH states on a torus with L\ = 5.3. 
The energy gaps of the finite systems with the Trugman- 
Kivelson-type potential are obtained by exact diagonal- 
ization, and extrapolation to the infinite-size system is 
done using the minimum-square method. The energy 
gaps of the corresponding spin chains are also calculated 
in a similar way. As shown in Fig. [51 we find that the 



energy gaps of the effective spin chains and the origi- 
nal FQH systems decrease with increasing p. We should 
note that energy gaps in the effective spin Hamiltonian 
do not always correspond with those in the original sys- 
tems, since the Hirbert space is limited. However, these 
gaps give upper bounds of the original gaps, at least. 
We also calculate the correlation function (S z (0)S z (n)) 
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of the total spins of the unit cell S'(n) = Efctp+i 
of the effective spin chains using LTEBD22 (see Fig. 
This numerical method enables us to calculate correla- 
tion functions in infinite-size systems without extrapola- 
tions. From Fig. 3, we can read off that the correlation 
functions vanish exponentially as functions of the dis- 
tances, regardless of whether the effective Hamiltonian 
includes the static terms or not. As shown in Table UU 
the correlation length of the effective spin model for the 
v = p/ (2p + 1) FQH state, £, where the length of a unit 
cell set to be unity, increases as p increases. These prop- 
erties of energy gaps and correlation functions are similar 
to those of the integer-,!? quantum spin chain where the 
Haldane gaps decrease and correlation lengths increase 
as S is increased. 



V. CONCLUSION AND DISCUSSION 

We have shown that FQH states in m = 1 Jain series, 
filling factor v = pj (2p+l) and v = p+l/(2p+l), beyond 
the thin-torus limit can be mapped to S = 1 effective 
spin chains with p sites in a unit cell. By numerically 
analyzing the energy gaps and the correlation lengths, 
we have shown that those effective spin chains and the 
original systems behave similarly in their p dependence. 
From these results, we point out that these effective spin 
chains have similar properties to those of S = p inte- 
ger Heisenberg chains with Haldane gaps which decease 
as S is increased. The above results give one of the ex- 
planations regarding the relationship of the hierarchy of 
FQH states and Haldane conjecture on the quantum spin 
chains. 

In the present spin- mapping for v = p/(2p+ 1), it was 
essential that the most relevant pair hopping process is 
given only by V21. Therefore, the Laughlin series v = 1/q 
with q > 5 cannot be treated in the same way, where 
contribution from the longer range hopping terms are 
important. Analysis for these states will be discussed 
elsewhere. 24 



VI. ACKNOWLEDGMENT 

We thank Emil J. Bergholtz for many helpful discus- 
sions. M. N. and Z.-Y. W. acknowledge support from 
the Global Center of Excellence Program "Nanoscience 



Filling factor v of the original system 


Correlation length £ 


1/3 


0.29 


2/5 


1.10 


3/7 


1.63 



TABLE II. Correlation length of the spin-spin corelation func- 
tion of the effective spin chains. The lengths of unit cells are 
set to be unity. 
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Appendix A: Proof of Eq. (flO|) 

Let us show that the all states in our model are gener- 
ated from the root state \^> r ) by applying only Uj defined 
by ijHJ) without using Uj, as written in (fTt)|). From Eq. ([9]), 
a configuration is supposed to be written as 



\*fa,...])=UtU r --\* r ), 



(Al) 



where the parts • • • are the products of U. The following 
conditions for operator Uj can be confirmed with simple 
calculations, 



Ulll,=0, (\k-j\<2), 
[Ul,Uj}=0, (\k-j\>3). 



(A2a) 
(A2b) 



Due to the space inversion symmetry, we only need to 
consider the nonvanishing case (k — j = 3). We get 



U-^U-j — c] +3 c] +6 c] +1 c] +2 c :)+5 c J - +4 c :(+3 c J - . (A3) 



Since ■ • • |* r ) ^ |- • • 111 ■ ■ •) in Eq. ([AT]) . (fj +3 Uj ■ ■ 
vanishes. Thus the above lemma has been proven. 



Appendix B: Proof of Eq. (fT5]) 



To prove the property described by Eq. (|15|) . we need 
the following lemma in a system with periodic or open 
boundary conditions: 



+1 ■ ■ 



■u k - 



■Uk - ■ ■ Uk+i 



0, (Bl) 



and 



where • • ■ means products of U. Since only Uj-3 
Uj+2p 3 +2 may reduce the number of particles in the un- 
derlined part of Eq. (??), we consider a case when one 
particle goes out from the unit cell by operating Uj-3. 
This situation is possible when a particle is located on 
the j-th site, 



!*> = !•• -oi. 



(B2) 



This state is generated after Uj-2 or Uj-i has been oper- 
ated to |\EV)- Since Uj-i increases the number of particles 
in the underlined part, we should only consider the case 



(B3) 



Now let us prove our lemma (jBlj) . We suppose N s to be 
the number of sites in the current system and \k—j\ < N s . 



7 



The operator Uj defined in (|8]) has the following proper- 
ties which can be confirmed with simple calculations: 



U j U k = U k Uj, (\k-j\^2), 
UjU k = 0, (4> \k-j\^2), 
UjUkUj = 0, 



(B4a) 
(B4b) 
(B4c) 



and 



UjUj+2 ■ ■ ■ Uj+2n 

=hj+2{nj+3 — 1) • • • flj + 2n(nj+2n+l 

X S'S'+l^j+2n+2^2+2r»+3' 
Uj+2n ' ' ' Uj+2Uj 

=nj+2n+l(nj+2n — 1) • ■ ' ?\j"+3 (%+2 

X S , +2n+3S , +2n+2S'+lS'- 



1) 



(B5a) 



1) 



(B5b) 



We find the following condition in the systems with 
periodic or open boundary conditions: 



Uk+i^JJk ••• = ••• U k -_jU k+1 • • • = 0, 



(B6) 



where the underlined parts do not include 
U k ,U k +i,U k -i, and U k +2- In other words, if we 
operate both U k and U k +i to |^ r ), then it vanishes. 
Let us prove this lemma. Because of the space inver- 
sion symmetry we only need to prove the condition 
• • • U k+ \-_^Jj k ■ ■ ■ = 0. When the system has open 
boundary conditions, wc just need to consider the 
following case: 



[U kl U l ]=Q(k^l±2)^ 

[fe,M = 0(Me^; |2s - 2* — 1| <N S 



(B7) 



where Eq. (|B4p has been used. It is obvious that the 
operators U k with even k and Ui with odd I commute 
each other, so that for the operators U in the underlined 
part, we can move all C/fc+ dd to the right of U k and all 



t4+cven to the left of Uk+i- Then we can change the 
order of U in the following way: 



' U k+ \-_j_JJ k 



■u k+x u k - 



0. 



(B8) 



Then it follows from Eq. (jB4bj) that (jBlj) has been proven 
for open boundary systems. 

For periodic boundary systems, we have to consider 
the relation U k — U k +\N S with A 6 Z. If N s is even, we 
can get Eq. (|B8I) in a similar way of the open boundary 
systems. On the other hand, if N s is odd (N s = 2M — 
1, M S N), we consider the following approach. First, we 
verify the following relation which is non-trivial in this 
case: 



[U k ,U k +2n] =0, (2<n<N s ). 



(B9) 



For 2 < n < (N s — l)/2, this relation can be obtained 
from Eq. ([B4|. For (iV s - l)/2 < n < N s , it follows from 
the relation U k +2n = U k+ 2(n-M)+i that the indices of U 
are always odd and less than N s . Therefore Eq. (|B9|) is 
satisfied from Eq. (IB4[) . Second, it follows from Eq. (|B5[) 
that the following relation is satisfied for n> M: 



U k +2n ■ ■ ■ U k+ 2U k ~ U k -2 



■ U k -2U k = o. 



(B10) 



Using Eqs. (|B9|) and (IBlOj) . if -_JJ k ■■■ is not 0, the op- 
erators U k -2 S ■ • • U k -iU k -2 in the underlined part can be 
moved to the left side of U k as 



■_JJ k ■■■ = 

{U k -2s ■ ■ ■ Iffc-4t4-2)t4(remaining U im 



(Bll) 



where < s < M. Finally, we consider the 
case • • • U k+ i-\ Ne ^JJ k ■ ■ ■ ^ '-'Uk+i-XNsUk'" = 
■ ■ - Uk+iUk - ■ ■ = 0. Comparing Eq. (|B11I) with 
• • • U k +i-xN 3 : • ~Uk • • ■ , the relation among s, A, and N s 
should be 



2s = AiV, - 1. 



(B12) 



Since iV s is odd, A is also odd. Therefore s > M, and 
Eq. (|B10|) yields • • • U k+1 -_J) k • • ■ = 0. Thus the lemma 
(|Bip has been proven for all cases. 

The lemma (jBlj) can also be easily proven by using 
the apagogical argument which assumes that one particle 
cannot move two sites in one direction from the initial 
condition. 
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